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ABSTRACT 

We propose an evolutionary path for prestellar cores on the radius-mass diagram, 
which is analogous to stellar evolutionary paths on the Hertzsprung-Russell Diagram. 
Using James Clerk Maxwell Telescope (JCMT) observations of L1688 in the Ophi- 
uchus star-forming complex, we analyse the HCO"'" (J=4— )-3) spectral line profiles of 
prestellar cores. We find that of the 58 cores observed, 14 show signs of infall in the 
form of a blue-asymmetric double-peaked line profile. These 14 cores all lie beyond the 
Jeans mass line for the region on a radius- mass plot. Furthermore another 10 cores 
showing tentative signs of infall, in their spectral line profile shapes, appear on or just 
over the Jeans mass line. We therefore propose the manner in which a prestellar core 
evolves across this diagram. We hypothesise that a core is formed in the low-mass, 
low-radius region of the plot. It then accretes quasi-statically, increasing in both mass 
and radius. When it crosses the limit of gravitational instability it begins to collapse, 
decreasing in radius, towards the region of the diagram where protostellar cores are 
seen. 
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> . 1 INTRODUCTION 

' Star formation in molecular clouds occurs within pre- 
?H ' protostellar (or prestellar for short) cores, which 
^ are Rra vitationally bound starless cores within the 



(|Ward- Thompson et al.l|l999l ). probed 1.3 mm continuum 
emission from nine prestellar cores with the IRAM tele- 
sco pe and studied their de n sity p r ofiles. Papers IV - VII 
(Je ssop fc Ward-ThompsonI l200ll : IWard- Thompson et all 



clouds jWard-Thompsoii et"all Il994 : Andre et a,1.1 



Ward- Thompson et al 
Ward- Thompson et al 



.19991 : 
I2OO7I ). 



Ward-Thompsoni " 



199e 



2002 
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Prestellar cores represent the initial conditions for pro- 
tostellar collapse. Their large-scale physical and chemical 
properties have be e n extensively studi e d. See, for example , 
iLoren et a l! (199d '):lMotte et al] ll 19981 ): I Andre ct al.l (|2007l ') 
and also DiFrancesco et al. 1 20071 . and references therein). 
However, the m anner in which they evolve i s still a matter 
of some debate (|Ward- Thompson et al.ir2007l . and references 
therein). In this paper we attempt to address the subject of 
prestellar core evolution and quantify this evolution with 
an empirical diagram analogous to the Hertzsprung-Russell 
Diagram. 

In this series of papers we have been investigat- 
ing the properties of isolated prestellar cores. The term 
'isolated' refers to the idea that these cores are not 
int eracting with each other (|Andre et all |2007| ). Paper 
I l|Ward- Thompson et al.l Il994l) marked th e discovery of 
prestellar cores. Paper II (|Andre et al.lll996l ') and Paper III 



SCUBA and Spitzer observations respectively to determine 
the density and temperature profiles of prestellar cores as 
well as estimating othe r properties such as t heir masses and 
hfetimes. Paper VIII (|Simpson et al.l I2OO8I ') used SCUBA 
data to measure the masses and radii of prestellar cores in 
L1688 to demonstrate that the stellar Initial Mass Function 
(IMF) has its origins in the C ore Mass Function (CMF). 
Paper IX (|Nutter et al.l |2009| ) studied prestellar cores in 
the far-infrared, using Akari observations, and demonstrated 
that such observations could be dominated by temperature 
gradients. In this paper we look at infall in prestellar cores. 

The Ophiuchu s star -forming region is located at a dis- 
tance of 139±6 PC llMama ick 2008) and is a site of low-mass 
star formation (|Wilking &: Ladal il983l.The region consists 
of two main clouds, L1688 and L1689, which have ex tended 
streamers leading out to distances of around 10 pc (|Lorenl 
Il989l ). Specifically, it is the more massive of the two clouds, 
L1688, that is studied in this paper, and that is generally 
known as the Oph main cloud. Very high star formation 
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Table 1. Table of regions showing date of observation, position (in Right Ascension, Declination), size (in arcseconds), position angle of 
observed HARP maps in this investigation. Right Ascension and Declination of map off positions also shown. All dates of observations 
are from 2009. 



Map 


Dates of 


R.A. 


Declination 


Width 


Height 


PA 


Off Pos. R.A. 


Off Pos. Declination 




Observations 


(2000) 


(2000) 


(") 


(") 


(°) 


(2000) 


(2000) 


Oph-A 


25-26 Apr, 01 May 


16''26™24.03= 


-24°23'43.0" 


690 


349 


45 


16''25™52.99= 


-24°25'23.0" 


Oph-AN 


12,20 Jun 


16''26™44.07" 


-24°18'53.7" 


420 


348 


50 


16''27™05.03" 


-24°24'58.3" 


Oph-B 


15,17,22 May, 02 Jun 


16'^27™22.93= 


-24°28'17.9" 


550 


348 


150 


16'^27™05.03= 


-24°24'58.3" 


Oph-C 


08,10-11 Jun 


16'^26™55.78= 


-24°34'16.5" 


420 


348 


50 


16''27™05.03= 


-24°24'58.3" 


Oph-E,F 


03,04,06-08 Jun 


16''27™15.76= 


-24°39'52.4" 


600 


348 


40 


16''27™07.82= 


-24°44'14.3" 



rates have been measured here wi th 14-40% o f the molecu- 
lar gas being converted into stars (|Vrballl977l ). 

The distance t o Oph iuchus is a matter of recent de- 
bate. iLoinard et al.l (|2008h use radio-emission parallax ob- 
servations from two young stars in Ophiuchus to c alculate a 

J irelim inar y distance measuremen t of 12014 2 P'-- lMamaie"5 
20081 ) and lLombardi et al.l (|2008l ) combine extinction maps 
with parallaxes from Hipparcos and Tycho. These studies 
find distances of 139±6 pc and 119=b 6 pc r espectively. For 
this analysis we adopt the iMamaiekl (|2008l ) value because 
it uses direct parallax measurements of a larger number of 
objects known to be associated with L1688. 

The Oph main cloud has been observed in 
many wavelengths from the visible to the submil- 



hmetre ((Montmcrlc et al.l Il983l: IWilkine et all 


198C; 


Loren et al.l 1990i: , Andre et al. 19921: Greene & Youne 


1992; 


Barsonv et al.l Il997l: iMotte et al.l 19981: 


Johnstone et al. 


20001: IJessop fc Ward- Thompson' '200l'; 


Johnstone et al. 


20041: [Stanke ct al. 2006: Enoch ct al. 2008: Jorgcnscn ct al. 


2008: IPadEett et all l2008l: iFriesen et al.l |2009|). Because 



this. 



its properties (e.g. mass, density, temperature, 
are very well known and it is therefore a good 



of 

distance 
plac e to study pres tcUar cores. 

ISimpson et al., ( 2008) used SCUBA archive data to pro- 
duce the deepest 850-^m map of L1688. This was used to 
determine the masses and radii of the prestellar cores in the 
cloud and to produce a core mass function (CMF). It was 
shown that the CMF can be fitted to a three-part power law 
consistent with the form of the stellar initial mass function 
(IMF). The form of the IMF was hypothesised to be a direct 
result of the CMF. A similar result was seen for t he Orion 
star-forming region bv lNutter fc Ward-ThompsonI (2007). 

In this paper we present HCO^ (J=4— >-3) observations 
of LI 688 taken using the Heterodyne Array Receiver Pro- 
gramme (HARP) instrument on the James Clerk Maxwell 
Telescope (JCMT). HCO+ (J=4^>3) was chosen because it 
is partially optically thick, which is a requirement for any 
line to show the characteristic profile for infall described in 
Sectioi JT^ In addition, HCO^ (J=4— >-3) can be observed at 
high spectral and spatial resolution with the JCMT. 



2 OBSERVATIONS 

HARP is a spectral-imaging receiver for the JCMT, op- 
erating at submilli metre wavelengths , in t he frequency 
range 325-375 GHz (|Smith et al.ll2003l . |2008| ). The HCO+ 
(J=4— >-3) line has a rest frequency of 356.7 GHz. HARP 
works in conjunction with ACSIS - Auto-Correlation Spec- 



tral Imaging System l|Dent et al.ll200Gl ') - a back-end corre- 
lator. It provides a spectral resolution of 0.026 kms~^, and 
a spatial resolution of 14 arcseconds. 

The observations here consist of 40 hours of data taken 
on multiple nights from April to June 2009. Five fields in 
L1688 were mapped, and these are outlined in Table [1] The 
maps were taken in raster position-switch mode. In this 
mode spectra are taken on the fiy, with the telescope scan- 
ning in a direction parallel to the sides of the map. At the 
end of each scan row, the array is displaced by one array- 
width perpendicular to the scan direction. Another map is 
created by scanning in a perpendicular direction to the first. 
These can be combined in a 'basket- weave' pattern to create 
a map with minimal striping due to differences in receptor 
noise. Separate fixed off-positions were used for each region. 
These are given in Table [T] These positions wer e selected as 
having no 850/im emission (jSimpson et al.ll2008l ). They were 
also checked for HCO^ (J=4— ^3) emission using a 60-second 
position-switched 'stare' observation. 

The data were reduced using t he StarlinkQ project soft- 
ware KAPPA JCurrie et al.! 20081, and references therein) 
and SMURF (jjenness et al.l l2008) routines. The data were 
regridded to a Nyquist-sampled map. It was found that in 
order to produce the optimal map from the input data, the 
worst quality data had to be discarded before the map was 
reconstructed. This was done by measuring the root-mean- 
square (RMS) noise value of each of the spectra in the input 
data and masking those with the highest RMS before gener- 
ating the map. An iterative method was used to determine 
the optimal RMS value to apply as a mask. In each step of 
the iteration, the data were masked using a different limit- 
ing RMS before being reconstructed into a map using the 
MAKECUBE routine (|jenness et al.ll2008l ). The RMS in the 
reconstructed map was then measured. The limiting RMS 
was varied to produce a reconstructed map with the lowest 
RMS. This map was taken to be the optimal map. 



3 RESULTS 
3.1 Maps 

Figure [T] shows the integrated intensity map towards L1688 
in HCO^ (J=4^3). A number of sources are seen in the 
map and the correlation between this integrated HCO"'" 
(J =4^3) map an d the s ubmillimetre continuum map shown 
in ISimpson et al.l (|2008l ) is consistent, showing that we are 
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Figure 1. Integrated intensity map of HCO+ (J=4-i>3) toward L1688. Contours at 1.0, 2.0, 3.0, 4.0 and 5.0 Kkms"!. The scale bar 
shows integrated intensity in Kkms~^. All temperatures refer to antenna temperature (T^). 



tracing the same material. The main regi ons of the cloud ar e 
labelled according to their usual names (jMotte at al.lll99^ . 

The RMS noise map for the HCO+ (J=4->-3) data is 
shown in Figure O Excellent weather conditions gave a very 
low RMS noise of less than 0.15 K in antenna temperature, 
T^, for maps of Oph-A, E and F. The noise in Oph-B in- 
creases to around 0.2 K toward the edges but most cores 
are still within the 0.15 K contours. Poorer weather dur- 
ing mapping of Oph-C resulted in a more mixed final map, 
with some portions mapped to 0.15 K and others to 0.2 K. 
The lower integration times on Oph-AN mean that the RMS 
noise level is higher in this region. 

Figure shows the same region as Figure [T] but with 
a narrower range of velocities. This figure has been split 
into red (4.0-5.0 kms~^), green (3.0-4.0 kms~^) and blue 
(2.0-3.0 kms~^) velocit y channels, centred around the sys- 
temic velocity of LI 688 (|Lorenll 19891 ). These colour channels 
were selected to encompass the velocity widths of the ob- 
served spectral line profiles. On this version of the map it is 
clearly seen that Oph-A is moving toward us relative to the 
rest of L1688. Both Oph-A and B show evidence for HCO+ 
(J=4— >-3) outflows and a wide range of velocities are seen 
in both regions. Oph-E and F appear much redder and are 
moving away from us relative to the other regions. Figure [4] 



shows the same data as Figure displayed as three separate 
images. 

Emission is dominated in the in tegrated intensity 
map by the bright Oph -A region (e.g. DiFrancesco et al.l 
|¥004), con taining SMI (|Wa.rd-Thompson et all Il989l ') and 
VLA1623 (|Andre et al.lll99^ ). respectively the prototypical 
prestellar core and Class object. Both of these sources are 
highlighted by tightly packed contours in Figure [T] A great 
deal of structure can be made out in Oph-A. An arch of ma- 
terial is seen traced out just northwest of the brightest part 
of the region - the northern part of which is moving toward 
us relative to the southern part. Just south of this arch, a 
line of sources make up a partial loop that curves around 
the western edge of the Oph-A region stretching through all 
the velocity channels. 

The Oph-B region (e.g. iFriesen et al.|[2005l ) is also very 
active and a distinct semi-circular loop can be seen on its 
south side. Both this loop and those in Oph-A have approx- 
imate radii of curvature of O.lpc and have kinematic veloci- 
ties of the order of 1-2 kms~^. The northern part of Oph-B 
is seen moving toward us relative to the rest of the region . 
This vel ocity gradient was also visible in ^^CO (iLorenll 19891 1 ■ 
DCO+ (|Loren et al.lll990l ). and N2H+ (|Andre et al.ll2007l ). 

Oph-C appears to be faint in the maps. The weak 
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Figure 2. HCO+ {J=4— >-3) noise map for L1688. Contours at 
0.15, 0.20 and 0.25 K. The scale bar shown is also in units of K. 
All temperatures refer to antenna temperature (T^). 



SCUBA sources measured here by ISimpson et al.l (|2008l ) 
seem to have little HCO"'' (J=4^3) associated with them. A 
large, ring-like structure is seen, in the red channel, making 
up the Oph-E region and this in turn leads into a filament 
that moves into the green channel as it joins the Oph-F re- 
gion. Only the northern-most part of the Oph-F region is 
seen in these maps. 



3.2 Spectral Line Profiles 

Prestellar cores undergoing collapse, will have excitation 
temperature profiles that are peaked toward their centre. 
In addition, we preferentially see material that is closer to 
us. Material on the near-side of a collapsing core is there- 
fore preferentially seen as less excited than material on the 
far side. This means that the relatively blue-shifted mate- 
rial is seen to the observer as having a higher excitation 
temperature. The result is that collapsing cores display an 
asymmetric spectral line profile, with the peak from mate- 
rial at the blue-shifted end being more prominent than the 
red-shifted peak. Such a profile shape is referred to as a 
blue-asymmetric double-peaked profile a nd can be used as 
a sign post that collapse is occurring (e.g. IZhoulll992l : lEvansI 
1 19991 ). 

We categorise all of the observed prestellar cores based 
on their HCO"'' (J=4^3) spectral line profile shapes: cores 
with blue-asymmetric double-peaked profiles; cores with 
faint or only possible blue-asymmetric double-peaked pro- 
files; cores with red-asymmetric double-peaked profiles; 
cores with single-peaked profiles and cores with unclear 
profile shapes. Observations of these cores using HCO"'" 
(J=4^3) are not sufficient to determine whether the core is 
undergoing collapse. Additional measurements of the cores, 
using an optically thin line tracer, are also required to de- 
termine that these double peaks are a spectral, rather than 
physical effect. For these data, we have included optically 
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Figure 9. HCO"'" (J=4^3) spectral line profiles for cores 
mapped for this study using HARP. Cores shown here have 
red-asymmetric do uble-peaked profiles. Systemic core velocities 
l|Andre et al.ll2007h are shown as dotted vertical lines when avail- 
able. 



thin N2H+ line centres from I Andre et all (|2007h . These data 
were available for the majority of the observed cores and are 
shown on the 3x3 pixel maps in the Appendix. 

Table [2] shows the cores mapped in this study, 
grouped into regions, along with their HCO^ (J=4^3) 
profile shapes. The masses and radii for these cores from 
ISimpson et al.l (|2008l ) are also shown. Figures [5] - [9] show 
the HCO^ (J=4— >-3) spectral line profiles for the same cores, 
grouped by profile shape. 



4 DISCUSSION 

If we assume that the prestellar cores are at a con- 
stant gas-kinetic temperature and are th ermally supported 
(|johnstone et al.l |2000| . I2OOII : iKirk et al.l 12007] '). then they 
will have a radius-mass relationship determined by the ex- 
ternal pressure on the m. The exact solutions ar e given by 
Bonnor-Ebert spheres (|Ebertlll955l : lBonnorlll956l '). This is in 
contrast to the cloud at large, where densities a re lower and 
turbulent motions are greater (|Kirk et al.ll2007t ). 

Assuming an approximately uniform and constant tem- 
perature for all the gas in the cloud allows us to calculate 
the Jeans Mass ( |Jeansiil902 ') . Below the Jeans Mass, the in- 
ternal forces in the core are dominated by thermal pressure. 
When above the Jeans mass, the core becomes unstable and 

collapses under its own self-gravity. 

Figure [TD] shows the cores detected by ISimpson et al.l 

(|2008t ) at 850/im with SCUBA, plotted in terms of 
radius versus mass. Here the radius is defined us- 
ing the 3(7 con tour of the signal-to-noi s e ma p of 
the region - see iNutter fc Ward-ThompsonI l|2007l ) and 



iNutter fc Ward-ThompsonI l|2007l )." The crosses mark the 
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Table 2. Table showing the Oph cores f rem f simpson et al mapped in this study by HARP. For each core, the radius and mass 

are given along with the spectral profile shape from our HARP data. Radii are rounded to the nearest 100 AU. Bl ue-asymmetric and 
red-asymmetric double-peaked profiles are listed as BAD and RAD respectively. Whether the core was included in the Motte et al.' l|l998l l 
— MAN98 — is shown for reference, along with an indication of whether the core was seen as a BAD profile in . Andre c t al., ( 2007. ) - A07. 



Core 


R.A. 


Declination 


Radius 


Mass 


Profile 


MAN98 


A07 


ID 


(2000) 


(2000) 


(AU) 


(Mq) 


Shape 


Core 


BAD 


A-MMl 


16*^26^23.00= 


-24°23'34.46" 


1000 


0.09 


Single 


Y 




A-MM2-3 


16*" 26™ 23. 80= 


-24°24'09.95" 


1200 


0.16 


BAD (Possible) 






A-MM4 


16'^26™23.03= 


-24°21'59.81" 


1000 


0.48 


BAD 


Y 




A-MM5 


16'^26™26.66= 


-24°22'28.60" 


1300 


0.82 


RAD 


Y 




A-MM6 


16'^26™28.21= 


-24°23'00.10" 


1500 


2.10 


RAD 


Y 




A-MM7 


16'"26™30.05= 


-24°22'17.93" 


1400 


0.92 


BAD 


Y 




A-MM8 


16'^26™31.80= 


-24°24'50.00" 


1000 


3.18 


Single 


Y 




A-MMll 


16'' 26™ 32. 74= 


-24°26'14.40" 


1600 


0.64 


Single 






A-MM16 


16''26™36.26= 


-24°28'12.84" 


1000 


0.03 


Single 






A-MM17 


16''26™34.77= 


-24°28'08.10" 


1000 


0.03 


Single 






A-MM18 


16''26™43.73= 


-24°17'25.74" 


1700 


0.87 


Single 






A-MM21 


16''26™31.64= 


-24°18'38.05" 


1000 


0.21 


Single 






A-MM22 


16''26™31.45= 


-24°18'52.05" 


1200 


0.11 


Unclear 






A-MM23 


16''26™07.89= 


-24°20'30.50" 


1600 


1.50 


BAD 






A-MM26 


16''26™ 15.40= 


-24°25'32.50" 


1200 


0.43 


BAD (Possible) 






A-MM27 


16''26™ 13.85= 


-24°25'25.16" 


1400 


0.34 


BAD (Possible) 






A-MM30 


16''26™09.63= 


-24°19'43.25" 


1500 


1.28 


Unclear 






A2-MM1 


16'^ 26™ 11.73= 


-24°24'54.16" 


1300 


0.33 


BAD (Possible) 


Y 




A3-MM1 


16''26™ 10.07= 


-24°23'11.00" 


1200 


0.33 


Single 


Y 




A-N 


16'' 26™ 22. 74= 


-24°20'00.00" 


1000 


0.21 


Single 


Y 




A-S 


16''26™42.69= 


-24°26'08.05" 


1000 


0.01 


Single 






SMI 


16''26™27.73= 


-24°23'58.17" 


1800 


7.35 


RAD 


Y 




SMIN 


16''26™27.93= 


-24°23'31.67" 


1100 


2.91 


Symmetric 


Y 




SM2 


16''26™29.41= 


-24°24'26.69" 


2200 


5.97 


RAD 


Y 


Y 


VLA1623 


16''26™26.74= 


-24°24'30.00" 


1600 


2.93 


RAD 






Bl-MMl 


16''27™09.32= 


-24° 27' 43.73" 


1100 


0.02 


Unclear 


Y 




B1-MM2 


16''27™12.14= 


-24°29'34.42" 


1200 


0.63 


BAD 


Y 




B1-MM3 


16''27™ 12.68= 


-24°29'38.67" 


1300 


2.61 


BAD 


Y 




B1-MM4 


16''27™ 15.32= 


-24°30'36.82" 


2000 


2.10 


Unclear 


Y 




B1-MM7 


16''27™18.72= 


-24°30'24.64" 


1000 


0.11 


Unclear 






B1B2-MM1 


16''27™ 12.44= 


-24°27'30.31" 


1800 


0.22 


Unclear 


Y 




B2-MM2 


16''27™ 19.97= 


-24°27'23.45" 


1400 


0.88 


BAD (Possible) 


Y 




B2-MM4 


16''27™ 24.50= 


-24° 27' 46.30" 


1000 


0.96 


BAD 


Y 




B2-MM5 


16''27™ 24.74= 


-24°27'29.29" 


1500 


1.14 


BAD 


Y 




B2-MM6 


16''27™25.27= 


-24°26'50.48" 


2300 


1.62 


Single 


Y 




B2-MM7 


16''27™26.66= 


-24°27'39.72" 


1700 


0.62 


BAD (Possible) 


Y 




B2-MM8 


16''27™27.96= 


-24°27'06.85" 


1300 


1.63 


Unclear 


Y 




B2-MM9 


16''27™28.74= 


-24°26'40.65" 


1400 


1.37 


BAD 


Y 




B2-MM13 


16''27™32.95= 


-24°26'03.16" 


1500 


0.49 


BAD (Possible) 


Y 




B2-MM14 


16''27™32.53= 


-24°26'27.78" 


1500 


1.26 


BAD 


Y 




B2-MM16 


16''27™34.49= 


-24°26'11.90" 


1200 


1.31 


BAD 


Y 


Y 


C-MM2 


16''26™59.10= 


-24°33'50.15" 


1000 


1.11 


BAD 


Y 




C-MM3 


16''26™58.80= 


-24°34'23.40" 


1900 


1.42 


Unclear 


Y 




C-MM5 


16''26™59.77= 


-24°34'13.87" 


1100 


1.24 


BAD 


Y 


Y 


C-MM6 


16''27™01.58= 


-24°34'44.62" 


2400 


1.05 


BAD (Possible) 


Y 


Y 


C-MM8 


16'' 26™ 49. 09= 
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Figure 3. False colour image of HCO+ (J=4-5.3) toward L1688 in red (4.0-5.0 kms-^), green (3.0-4.0 kms^^) and blue (2.0-3.0 kms^^) 
velocity channels. 



prestellar cores. The cores with associated infrare d sources 
ident ified by Spitzer, as part of the c2d survey ( Enoch et al.l 
I2OO8I ). are deemed to be protostellar and are marked with 
asterisk s. We assume that these objects are C lass or I pro- 
tostars l| Andre et al.ll 19931 : lEnoch et al]|2009l ) , meanine that 
they have accreted only part of their final mass. In order to 
account for the protostellar mass not detected by SCUBA 



we use a mean IMF mass of 0.6Mq (|Kroupall2002l ) and thus 
increase the masses of these objects by 0.3Mq. This would 
be the correction needed if all of the objects lay exactly on 
the Class O/I borderline, defined as the point at which a 
protostar has accr eted half of its final main sequence mass 
(j Andre et al.ll 19931 ) . and hence is indicative of the correction 
needed. 



ited evolutionary diagram for prestellar cor 



16:27:00 30 
R.A, (2000) 
(a) 4.0-5.0 kms-l 



30 16:27:00 30 

R.A. (2000) 
(c) 2.0-3.0 kms-l 



000) 

kms-i 



I 



Figure 4. (a) Integrated intensity maps of HCO+(4-3) toward L1688 in red (4.0-5.0 kms~^), green (3.0-4.0 kms~^) and blue (2.0-3.0 
kms~^) velocity channels. The scale bar shows integrated intensity in Kkms~^. All temperatures refer to antenna temperature (T^). 



We propose the radius versus mass diagram shown in 
Figure fTOl to be the equivalent of an HR diagram for prestel- 
lar cores. Cores evolve by accreting mass quasi-statically and 
maintaining Bonnor-Ebert equilibrium. Upon reaching the 
Jeans instability, cores collapse to form a protostar, moving 
downward on this diagram. 

Protostellar cores lie beyond the Jeans limit, to the right 
on the graph, near the resolution limit of the observations. 
These cores have already accreted material and begun to 
collapse. Again these cores have been given an additional 
0.3Mq to account for the undetected mass. 

The hypothesis that this diagram illustrates the evo- 
lution of prestellar cores can be tested by looking at the 



position on this diagram occupied by collapsing prestellar 
cores. If this evolutionary picture is correct, then the blue- 
asymmetric double-peaked profiles should all lie close to, or 
beyond the Jeans mass. 

Table [5] lists the masses and radii of those starless 
cores detected by SCUBA, that were also detected in HCO'*' 
(J=4— >-3) by HARP. Their profile shapes and positions are 
also given. Figure [TT] is the same as Figure [TO] except that 
it shows only the cores listed in Tabled This time they are 
plotted according to their spectral line profile shapes. Those 
cores with blue-asymmetric profiles, indicating collapse, are 
marked with dark boxes. Cores showing possible infall char- 
acteristics are marked with hollow boxes. All other cores are 



8 Simpson et al. 




2 3 4 5 

velocity {kms ') 



velocity (Icms ) 



2 3 4 

velocity (kjTis ') 



2 3 4 

velocity (kms ') 



Figure 5. HCO"'' (J=4— >3) spectral line profiles for cores mapped for th is study usiiig HA RP. Cores shown here have characteristic 
blue-shifted, asymmetric double-peaked profiles. Systemic core velocities llAndre et al.ll2007f) are showrn as dotted vertical lines when 
available. 
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Figure 6. HCO+ (J=4— f3) spectral line profiles for cores mapped for thi s study using HAR P. Cores shown here display profiles with 
possible blue-shifted, asymmetric double-peaks. Systemic core velocities llAndre et al.|[2007l ) are shown as dotted vertical lines when 
available. 
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Figure 7. HCO"'" (J=4— >3) spe ctral line profiles f or cores mapped for this study using HARP. Cores shown here have single-peaked 
profiles. Systemic core velocities jAndre et al.|[2007l) are shown as dotted vertical lines when available. 



2.0 
1.5 
1.0 
0.5 
0.0 

2.0 
1.5 
1.0 
0.5 
0.0 

2.0 
1.5 
1.0 
0.5 



Bl-MMi 




I I I- 




1 1 r 



I I I- 



Bi-llM7(3'x) ' 




:-lv[Mf2(3'x) ' 



c-ivImioo'x) ' ' 




Ir 



I I Ml 



C-Iv1m8(3x') ' 




2 3 4 5 

velocity (kms ') 



2 3 4 5 

velocity (kms ') 



2 3 4 

velocity (kms ') 



2 3 4 5 

velocity (kms ') 



Figure 8. HCO+ (J=4— >3) spectral line profiles for cores mapped for this study using HARP . Cores shown here are those displaying 
profiles of no particular shape or low signal-to-noise. Systemic core velocities jAndre et al.ll2007l ) are shown as dotted vertical lines when 
available. 
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Figure 10. Starless and protostellar cores from lSimpson et al. I ||2008|'I plotted in terms of radius against mass. Starless cores are indicated 
with crosses. Protostellar cores (i.e. cores with infrared associations) are marked with asterisks. For error bars, see Figure [TT] See text 
for more details. 



marked as crosses. These cores may simply have the wrong 
physical conditions to produce a detectable blue-asymmetric 
profile in HCO^ (J=4— >-3). Figure [9] shows those cores that 
display red-asymmetric profiles. These are hkely caused by 
bipolar outfiows confusing the signature, especially in the 
case of VLA1623, which might also interfere with the neigh- 
bouring SMI and SM2. It is noted that this explanation does 
not necessarily account for the profile of the core C-N, which 
displays a very thin RAD profile. 

In Figures [TO] and [TT] the two dashed lines marked 
'Low Pressure' and 'High Pressure' show the Bonnor- 
Ebert relations for external pressu res of P/fc = 3 x 10^ 
and 12 x lO'' K cm~ ^ respectively (|Johnstone et all [2OO0I : 
IJorgensen et al.l l2008h . The solid line in Figure [11] shows a 
possible evolutionary track, indicating the general direction 
a core would move through the diagram. This track initially 
follows a Bonnor-Ebert relation for an external pressure of 
P/fc = 6 X 10^ K cm~'^. We assume a temperature for the 
cores in Ophiuchus of 10 K (jStamatellos et aDl2007h . 

The hypothesis that a radius-mass plot can serve as an 
evolutionary diagram for prestellar cores is supported by 
these results. In Figure (TT] cores showing only weak signs of 
infall (hollow boxes) appear close to, or beyond, the Jeans 
mass line compared with those cores exhibiting characteris- 
tic infall signatures (solid boxes) which all appear beyond it. 



40% of the cores beyond the Jeans mass line show signs of 
infall; 60% if the possible BAD peaks are included. This sup- 
ports the conjecture that cores follow an evolutionary track 
in which they first move diagonally up the diagram, as they 
accrete material, and then move down it as they cross the 
Jeans mass limit and collapse. 

Using this model, for any given prestellar core on this 
diagram, an approximate protostellar mass could be derived 
by following a Bonnor-Ebert track diagonally upward to the 
Jeans mass line. From that point in the evolution, an effi- 
ciency for the collapse of the core needs to be assumed to 
reach an estimate of the final protostellar mass. If a core 
continues to accrete whilst collapsing, the gradient of its 
evolutionary track, after crossing the Jeans mass, will be 
shallower. 



5 CONCLUSIONS 

We propose that a radius-mass plot can serve as an evo- 
lutionary diagram for prestellar cores. Of the 58 prestellar 
cores observed with HARP, 14 showed signs of infall in the 
form of a blue-asymmetric double-peaked line profile. These 
14 cores all lie beyond the Jeans mass line for the region on 
such a radius-mass plot. Furthermore those cores showing 
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Figure 11. Same as Figure [TOl excluding tlie protostellar cores. Cores sliowing definite blue-asymmetric double (BAD) peaks in their 
spectral line profiles are shown as dark squares; cores with a possible blue-asymmetric double peak are shown as hollow squares; all other 
cores are shown as crosses. Systematic errors are shown in the top-left hand of the plot for distance to Ophiuchus, D, temperature of the 
cores, T and dust mass opacity, k. An error of 10% is estimated for the measurement of the radii of cores. This produces a corresponding 
error in mass. These error bars are shown on the plot. We propose this as an evolutionary diagram for prestellar cores, and indicate the 
manner of this evolution with an arrow. 



tentative signs of infall, in their spectral line profile shapes, 
appear close to, or beyond, the Jeans mass line. 

Placing prestellar cores in the context of this evolution- 
ary diagram shows how they can evolve into protostars in- 
side molecular clouds. This analysis could be performed us- 
ing data for other star-forming clouds, where both submil- 
limetre maps and spectral profile data are available. 

Future studies using instruments such as SCUBA-2 and 
Herschel may be able to probe further to the left on dia- 
gram, and see how larger, lower-mass objects - such as CO 
cores - relate to the gravitationally bound prestellar cores 
already seen in the continuum. Such objects appear to fol- 
low different radius-ma ss relations to those explored in this 
st udy - see Figure 4 of [Motte et al.l (|200ll ) and discussions 
by I Johnstone et al.l (|200(]| . I2001I ')" 
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APPENDIX A: 3X3 PIXEL HARP MAPS. 

Figures |XT] - 1X581 show the HCO+ (J=4^>3) spectral Une 
profiles for the same cores seen in Figures [5] -[9] These plots 
show the central 3x3 HARP pixels for each core. These plots 
also show the N2H"'" line centres from Andre ct al. (2007) as 
vertical dashed lines. These data are not available for all 
cores. 
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Figure A2. Map of HCO+ (J=4-^3) spectra in tlie core A-MM2-3. The central grid of 3x3 HARP pixels is shown here. 
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Figure A4. Map of HCO+ (J=4— ^3) spectra in the core A-MM5. The central grid of 3x3 HARP pixels is shown here. 
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Figure A5. Map of HCO+ (J=4— >-3) spectra in tlie core A-MM6. Tlie central grid of 3x3 HARP pixels is shown here. 
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Figure A6. Map of HCO+ (J=4-»-3) spectra in the core A-MM7. The central grid of 3x3 HARP pixels is shown here. 
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Figure A7. Map of HCO+ (J=4— >-3) spectra in tlie core A-MM8. Tlie central grid of 3x3 HARP pixels is shown here. 




Figure A8. Map of HCO+ (J=4->3) spectra in the core A-MMll. The central grid of 3x3 HARP pixels is shown here. 
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Figure AlO. Map of HCO+ (J=4-^3) spectra in the core A-MM17. The central grid of 3x3 HARP pixels is shown here. 
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Figure All. Map of HCO+ (J=4^3) spectra in tlie core A-MM18. The central grid of 3x3 HARP pixels is shown here. 
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Figure A12. Map of HCO+ (J=4-^3) spectra in the core A-MM21. The central grid of 3x3 HARP pixels is shown here. 
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Figure A14. Map of HCO+ (J=4-^3) spectra in the core A-MM23. The central grid of 3x3 HARP pixels is shown here. 
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Figure A16. Map of HCO+ (J=4-^3) spectra in the core A-MM27. The central grid of 3x3 HARP pixels is shown here. 
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Figure A17. Map of HCO+ (J=4^3) spectra in tlie core A-MM30. Tlie central grid of 3x3 HARP pixels is shown here. 
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Figure A18. Map of HCO+ (J=4— ^3) spectra in the core A-N. The central grid of 3x3 HARP pixels is shown here. 



A suggested evolutionary diagram for prestellar cores 23 



1 1 

rr 1 




1 


A 

1 


I 1 1 


1 

^^^^ 


A 




1 1 
^ 1 1 




1 


1 


1 1 1 




1 




1 1 




1 


1 
1 


1 1 1 


1 


1 
1 





2345 2345 2345 

velocity {kms'b velocity (kms'^) velocity (kms'b 



Figure A19. Map of HCO+ (J=4— >3) spectra in ttie core A-S. The central grid of 3x3 HARP pixels is shown here. 
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Figure A20. Map of HCO+ (J=4-^3) spectra in the core A2-MM1. The central grid of 3x3 HARP pixels is shown here. 
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Figure A21. Map of HCO+ (J=4^3) spectra in tlie core A3-MM1. Tlie central grid of 3x3 HARP pixels is shown here. 
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Figure A22. Map of HCO+ (J=4-^3) spectra in the core SMI. The central grid of 3x3 HARP pixels is shown here. 
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Figure A23. Map of HCO+ (J=4^3) spectra in tiie core SMIN. Tlie central grid of 3x3 HARP pixels is shown here. 
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Figure A24. Map of HCO+ (J=4-^3) spectra in the core SM2. The central grid of 3x3 HARP pixels is shown here. 
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Figure A25. Map of HCO+ (J=4^3) spectra in the core VLA1623. Tlie central grid of 3x3 HARP pixels is sliown liere. 

I 1 1 n 1 1 1 1 n 1 1 1 1 n 1 1 




J I U I I I I li I I I I U L 

2345 2345 2345 

velocity (kms"^) velocity (kms'^) velocity {kms'^) 



Figure A26. Map of HCO+ (J=4-»-3) spectra in the core Bl-MMl. The central grid of 3x3 HARP pixels is shown here. 
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Figure A28. Map of HCO+ (J=4-»-3) spectra in tiie core B1-MM3. The central grid of 3x3 HARP pixels is shown here. 
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Figure A30. Map of HCO+ (J=4-»-3) spectra in tiie core B1-MM7. The central grid of 3x3 HARP pixels is shown here. 
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Figure A32. Map of HCO+ (J=4-»-3) spectra in tiie core B2-MM2. The central grid of 3x3 HARP pixels is shown here. 
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Figure A34. Map of HCO+ (J=4-»-3) spectra in tiie core B2-MM5. The central grid of 3x3 HARP pixels is shown here. 
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Figure A35. Map of HCO+ (J=4->-3) spectra in tlie core B2-MM6. The central grid of 3x3 HARP pixels is shown here. 




Figure A36. Map of HCO+ (J=4-^3) spectra in the core B2-MM7. The central grid of 3x3 HARP pixels is shown here. 
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Figure A38. Map of HCO+ (J=4-»-3) spectra in tiie core B2-MM9. The central grid of 3x3 HARP pixels is shown here. 
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Figure A40. Map of HCO+ (J=4-)-3) spectra in the core B2-MM14. The central grid of 3x3 HARP pixels is shown here. 
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Figure A41. Map of HCO+ (J=4->-3) spectra in the core B2-MM16. The central grid of 3x3 HARP pixels is shown here. 
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Figure A42. Map of HCO+ (J=4-^3) spectra in the core C-MM2. The central grid of 3x3 HARP pixels is shown here. 
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Figure A43. Map of HCO+ (J=4^3) spectra in the core C-MM3. Tlie central grid of 3x3 HARP pixels is sliown liere. 
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Figure A44. Map of HCO+ (J=4-^3) spectra in the core C-MM5. The central grid of 3x3 HARP pixels is shown here. 
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Figure A46. Map of HCO+ (J=4-^3) spectra in the core C-MM8. The central grid of 3x3 HARP pixels is shown here. 
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Figure A48. Map of HCO+ (J=4-»-3) spectra in the core C-MMIO. The central grid of 3x3 HARP pixels is shown here. 
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Figure A49. Map of HCO+ (J=4->-3) spectra in the core C-MM12. The central grid of 3x3 HARP pixels is shown here. 
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Figure A50. Map of HCO+ (J=4->3) spectra in the core C-N. The central grid of 3x3 HARP pixels is shown here. 
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Figure A52. Map of HCO+ (J=4-^3) spectra in tlie core E-MM2b. The central grid of 3x3 HARP pixels is shown here. 
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Figure A53. Map of HCO+ (J=4->3) spectra in ttie core E-MM2d. Tlie central grid of 3x3 HARP pixels is sliown liere. 
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Figure A54. Map of HCO+ (J=4-^3) spectra in the core E-MM4. The central grid of 3x3 HARP pixels is shown here. 
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Figure A56. Map of HCO+ (J=4-^3) spectra in the core E-MM8. The central grid of 3x3 HARP pixels is shown here. 
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Figure A57. Map of HCO+ (J=4— >-3) spectra in tlie core F-MMl. Tlie central grid of 3x3 HARP pixels is shown liere. 
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Figure A58. Map of HCO+ (J=4— >3) spectra in the core F-MM2a. The central grid of 3x3 HARP pixels is shown here. 



